Quantum phase transition of the one-dimensional transverse field compass model 



O 

o 

> 

o 

(N 



Ke-Wei Sun^'^ and Qing-Hu Chen^-i^* 
^ Department of Physics, Zhejiang University, Hangzhou 310027, P. R. China 
^ Center for Statistical and Theoretical Condensed Matter Physics, 
Zhejiang Normal University, Jinhua 321004, P- R- China 
^ Institute of Materials Physics, Hangzhou Dianzi University, Hangzhou 310018, China 

(Dated: November 21, 2009) 

The quantum phase transition (QPT) of the one-dimensional (ID) quantum compass model in a 
transverse magnetic field is studied in this paper. An exact solution is obtained by using an extended 
Jordan and Wigner transformation to the pseudo-spin operators. The fidelity susceptibility, the 
concurrence, the block-block entanglement entropy, and the pseudo-spin correlation functions are 
calculated with antiperiodic boundary conditions. The QPT driven by the transverse field only 
emerges at zero field and is of the second-order. Several critical exponents obtained by finite size 
scaling analysis are the same as those in the ID transverse field Ising model, suggesting the same 
universality class. A logarithmic divergence of the entanglement entropy of a block at the quantum 
critical point is also observed. From the calculated coefficient connected to the central charge of the 
conformal field theory, it is suggested that the block entanglement depends crucially on the detailed 
topological structure of a system. 

PACS numbers: 05.70.Fh, 75.40.Cx, 73.43.Nq, 75.10.-b 
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I. INTRODUCTION 



The quantum compass model has been studied exten- 
sively in recent years due to the possible long range or- 
bital order and the quantum phase transitions (QPTs) [ll, 
0,11,0,1,11, 0,3. First, the model could be used to de- 
scribe the Mott insulators with orbit degeneracies. It de- 
pends on the lattice geometry, and belongs to the low en- 
ergy Hamiltonian originated from the magnetic interac- 
tions in Mott-Hubbard systems with the strong spin-orbit 
coupling[9i, jy]j. For simplicity, the ID quantum compass 
model is regarded as the coupling along one of bonds 
which shows an Ising type, but different spin components 
are active along other bond directions. It is exactly the 
same as the ID reduced Kitave model[lll, [H, [H, [II- 
The symmetry of the pseudo-spin Hamiltonian is much 
lower than SU(2). It is shown in the numerical results 
that the eigenstates are at least twofold degenerate or 
highly degenerate [l|,@. Recently, because of degeneracy 
in the ground-state (GS), the protected qubit can be im- 
plemented, a scalable and error-free scheme of the quan- 
tum computation can be designed in this simple model [1]. 

To shed some insights into this model, a exact solu- 
tion is clearly desirable. By applying Jordan and Wigner 
transformation to the pseudo-spin operators, Brzezicki et 
al. were able to map it into a spinless fermion model and 
determine the spectrum exactly [5]. More recently, the 
exact solution of ID period-two compass model has also 
been obtained by the present authors and a collaborator 
with a slightly different method^. In order to be useful 
for quantum information, the GS must be protected from 
local perturbations, but the spectrum for h — Q in gap- 
less in the thermodynamic limit. The extension to finite 
fields is a crucial step in the search for systems supporting 
naturally robust quantum information [Tsl . To the best 
of our knowledge, the ID compass model in the trans- 



verse magnetic field has not been studied so far, which 
may be also of fundamental significance. The symme- 
try of the system is further broken when the transverse 
magnetic field is applied. The behavior of the energy gap 
may be changed around the critical point in the thermo- 
dynamical limit and the degeneracy in the GS may be 
lifted, therefore the nature of the QPT may be altered 
in the presence of the transverse magnetic field. In order 
to address these questions, an exact solution to the field 
version is also clearly called for. 

Due to the recent progress in quantum information 
science, some concepts in quantum information theory, 
such as the fidelity, the fidelity susceptibility (FS), and 
the quantum entanglement have been extensively used to 
identify the QPTs in various many-body systems from 
the persMctive of the GS wave functions [H, [13, [H, 
[H, [il, [in, [12, [H, [H- Recently, it is proposed that 
the fidelity approach is a valuable tool to investigate 
novel phases lacking a clear characterization in terms 
of local order parameters (25i [26l |. With these effective 
tools and the finite-size scaling analysis of the FS, one 
can identify the universality class of the QPT in var- 
ious models [27l. [28j . Quantum entanglement is one of 
the most striking consequences of quantum correlation in 
many-body systems, and is recognized to be resource that 
enables quantum computing and communication[29| . It 
has shown a deep relation with the QPTjsl, l3l|. The 
entangled degree between any two nearest-neighbor par- 
ticles keeps the same for the translational symmetry, 
and its derivative may play the role of an order pa- 
rameter to characterize QPT at the critical point. In 
the context of QPTs, the quantum entanglement has 
been the subject of considerable interests in the various 
models m, HJ, [13, [H, [H . 

In this paper, we study the ID compass model in a 
transverse magnetic field for the first time. The exact 
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FIG. 1: (Color online) The odd {l,n;2,n} and the even 
{2,n;l,n+l} bonds are denoted in the regions of black dashed 
rectangle and blue dotted ellipse, respectively. 



solutions are obtained by using the method of mapping 
into a case with plural spin sites [s^. The GS fidelity, the 
FS, the concurrence and the block-block entanglement 
entropy are calculated. The behaviors of the spin corre- 
lations function are also given. The paper is organized 
as follows: In Section II, we describe the model and the 
scheme to obtain the exact solution in detail. The cal- 
culations of the fidelity, the concurrence and the block- 
block entanglement entropy are presented in Section III, 
where the scaling analysis is also performed. The correla- 
tion functions are analyzed in Section IV. The conclusion 
is given in the last section. 



boundary conditions (PBC) for pseudospins is employed, 
i.e. CTi^jv'+i = CTi^i. Note that the ID compass model 
without the magnetic field is just a special case of the 
alternating XY model 36] . 

In order to diagonalize the Hamiltonian (2), we use 
the extension of the Jordan and Wigner transformation 
for the case with plural spin sites [35]. An up-spin state 
is transformed to a one-fermion state, and a down-spin 
state to a zero-fermion state. The explicit mapping be- 
tween spin operators and fermionic operators are given 

by 



'{a\ + ai,„)(4,n - «2,n) 



- oJ 



(3) 



Here we denote the fermion creation operator with site 
number s and cell number n by a\ „. Then the Hamilto- 
nian (2) is transformed into the following form 



N' 



N' 



-a2.„a 



II. MODEL HAMILTONIAN AND EXACT 
SOLUTION 

The ID XX- YY model in a transverse magnetic field 
can be regarded as the structure of two pseudo-spin sites 
inside a unit cell. The Hamiltonian is given by[5| 

TV' N' 
n— 1 n— 1 

- J/3E <n4« - 2 T.^<n + ^2,„), (1) 
n— 1 n— 1 

where N — 2N' is the total number of the sites. Fig. 1 
shows the structure of interactions in Eq. (1). 

For /3 = 1, it becomes the ID compass model in a 
transverse magnetic field 



N' 



N' N' 

71—1 n—1 
, N' 



y 

2,n 



-9 



(2) 



where J denotes the strength of the nearest-neighbor in- 



teraction, a is the coupling parameter, a. 



are the 



Pauli matrix on cell n with site s = 1, 2, = 2A^' is the 
total number of the sites, and h is the applied magnetic 
field in the z direction. For convenience, the number 
of pseudospins N' is chosen to be even, and a periodic 



n=l 



4,na2,n) 



hN'. 



(4) 



The Fourier transformation of the fermion operators 
gives Qs^n = Z]p For convenience, 

the antiperiodic boundary condition (ABC) ai_Ar/-|-i = 
— oi^i is employed for the fermion operators. After these 
transformations, the new Hamiltonian H' now reads 

p 

-e'Pa2{p)a\{p)~e^Pa2{p)a,{-p)] 

-Ja'^[-a\{p)al{-p) + a\{p)a2ip) 
p 

-ai{p)al{p) - ai{p)a2{-p)] 

~hY,[4ip)ai{p) + 4ip)a2ip)] + hN', (5) 
p 

where p is the wave number in ABC which takes such 
values as p — ±jTr/N' , {j = 1,3,...N' — 1). The oper- 
ators al{p) and as{p) are the creation and annihilation 
operators of the fermion with site numbers s and wave 
number p, which satisfy the following anticommutation 
relations 

{a sip), alii)} = Ss,t6p,q, 
K(p),at(g)} = 0,{alip),aUq)} = 0. (6) 
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Then we find that the Hamiltonian H' is the sum of the 
foUowing independent operators W{p'): 

W{p') = -J[e-^^'a\{p')a\{~p') + e~^^'al{p')a^{p') 
-e'P a2{p')a\{p') - e'P' a2{p')ai{-p') 
+e^p' aU-p')a\{p') + e^^' al{^p')a,{~p') 



a2{-p')a\{-p') - e-^P a2{-p')a,{p')] 
-Ja[-a\{p')al{-~p) + a\{p)a2{p') 
~ai{p')al{p') - ai{p')a2{~p') 
-a\{-p')al{p') + a\{-p')a2{-p') 
-ai{-p')al{-p') - ai{-p')a2{p')] 
-h[a\{p')a,{p') + al{p')a2ip') 
+a\{-p')ai{-p') + al{-p')a2{-p')], (7) 

where p' = jTr/TV', (j = 1,3, ...TV' - 1). Note that 
[W{p'), W{q')] — 0, so we can solve the Hamihonian (5) 
in the space of p'. 

The parity in the Hilbert space of W{p') is conserved, 
therefore subspace with the even parity can be easily 
constructed in terms of the following 8 basis vectors 
|0),al(p')al(-p')|0),4(p')4(-p')|0),4(p')al(-y)|0), 
4(P')«1(-P')|0),al(p')4(p')|0),al(-P')4(-P')|0),and 
al{p')a\{p')al{-p')a\{-p')\0). While the 

subspace with the odd parity is obtained 
by combining the following 8 basis vectors 
4(p')|0),4(y)|0),4(-y)|0),4(-p')|0),al(-p')4b') 

4(y)lo),4(-y)4(p')4(p')|o),4(y)4(-p')4(-p')|o), 

and al{p')al{—p')a\{—p')\0). The parity of subspaces 
determines the boundary conditions. Indeed, the Bogoli- 
ubov vacuum has even (odd) numbers of a quasiparticles 
for ABC (PBC)[i]. For ABC, nonzero elements of 
the 8x8 Hermit matrix (even parity) for the reduced 
Hamiltonian W{p') are 

Wi,» = ~2h fori < i < 7, Wss = -4:h, 
W2,3 = W4,5 = ~Je-'P' - Ja, 
Wi^3 = Wi^s = Je^'P' + Ja, 
W2,A = VK3,5 = VKi,4 = W3,8 = -^e'^' - Ja. (8) 
The eigenvalues for this matrix are then easily derived 
A(i.2) (p') = _2h ± 2J[{h/jf + 1 + 2a cos(p') + a^Y'^ 
A(3,4)(p') = -2/i±2J[l + 2acos(p')+a2]i/2 
A(5-8)(p') = -2h. (9) 
The spectral functions e(p') are readily obtained 
g(i,2)(y) ^ ±2J[{h/jf + l + 2acos{p')+a^Y''^ 
g(3,4)(y) = ±2J[1 + 2a cos(p') (10) 
Actually, the Hamiltonian (5) can be decomposed as 



p' 



(11) 



where wjf'^ = e^''Kp')Tll\''\li^ (s' = 1,2,3,4) with r^^J^ 
the operator of fermionic quasiparticles , and the corre- 
sponding eigenvectors are ip^'^ \p') = 'yp'" ^|0). Then the 
GS energy and wave function are given by 

Eg = - ^2J[(VJ)' + l + 2acos(y)+a2]i/2, (12) 
p' 



p' 



(13) 



It should be stressed here that Eqs. (12) and (13) are 
valid for any value of h. Note also that the exact spec- 
trum is the same as that obtained by Brzezicki et al. for 
h — using a different method^. 

The energy gap can be readily obtained as A = 
2J[{h/jf + 1 ± 2a + a2]i/2 _ 2J\a ± 1|. It does not 
disappear in the presence of the transverse field even in 
the thermodynamic limit. The QPT driven by the trans- 
verse field will occur at (a = ±l,ft, = 0), which shows 
the second-order nature, similar to the QPT driven by 
interaction parameters [8[. 

For completeness, we will also briefly discuss the spec- 
tra based on PBC. For PBC, we need to solve the Hamil- 
tonian (7) in the odd numbers of a quasiparticles sub- 
space. The spectral functions are given by 



e(i^2)(y) ^ -J[(VJ)' + l + 2acos(p') 
-J[l + 2acos(p')+a^]^/^ 
e(3:4)(p') = -J[(/i/J)' + l + 2acos(p') + a']'/' 

e(5^6)(p') 

e(^-«)(p') 



-J[l + 2acos(p') 



v2ll/2 



J[{h/Jf + 1 + 2acos(p') + a2]i/2 
-J[l + 2acos(p') + a2]i/2 

J[{h/Jf + 1 + 2acos(p') + a^]^'"^ 



-J[l + 2acos(p') 



v2ll/2 



(14) 



Note that p' — Q and p' = tt in ABC must be treated 
separately and carefully. It is helpful to write down ex- 
plicitly the spectra for 4 pseudospin sites in the real 
space 

c^(i'2) = ±2J[[h/jf + l + a'Y'^ 



,(3,4) _ 



±2J[1 



.211/2 



.(5-8) ^ 
,(9,10) _ 



,2]l/2 
v2ll/2 



= Ja + J ±J[{h/JY + l-2c 

("'12) ^ Ja-J±J[{h/jf + l + 2c 

(13^14) ^ _-Ja_j±jp/j)2 + i_2a + a2]i/2 

(15^16) = -Ja + J±J[{h/jf + l + 2a + a^Y/^, 



(15) 



which include the spectra for both ABC [w^^ and 
PBC The GS energy for PBC can be written 
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FIG. 2: (Color online) GS energy as the function of the pa- 
rameter aioT h = 0.8, J = 1.0 and = 8 with PBC. 



where s'=9-16, and p' = 2j7r/iV', (j = 1, 2, ...^ - 1) 

The vaUdity of Eq. (16) is also confirmed by compar- 
ing with the direct numerical diagonalization of = 8 
pseudospin sites in real space. In Fig. 2, we present the 
numerical GS energy from PBC, and the results from Eq. 
(16) as well. It is clear that the present analytical results 
for the GS energy are in excellent agreement with the 
numerical ones. 

In order to show the correctness of the present method, 
we extend to study the celebrated ID Ising model in a 
transverse magnetic field, which Hamiltonian reads 

N' N' 

n— 1 n— 1 

, N' 



+ 2 IIKn+4n)• 
r^=l 

With ABC, the exact GS energy is derived as 
Eg = - Y.ih' + AJ^ ± 4J/1C0S ^f'\ 



(17) 



(18) 



where p = jn/N' , (j ~ 1,3, ...N' — 1), and the number 
of total sites is iV = 2N'. For PBC, the GS energy is 
written as 



G 



(/i' + 4j2±4J/icos^)i/2 



Py^Oom 



(19) 



= -2J - (4J2 + /j2)l/2 (2J > /i > 



where 



0), and p = 2j7r/A^', (j 



FIG. 3; (Color online) Fidelity Fmin = min[F{h,h + 
Sh),F{a,a + 5a)] in the a- /i plane for N' = 100 and J = 1.0 
with ABC. The second-order QPT points are obviously found 
at (a = -1.0, /i = 0) and (a = 1.0, h = 0). 



recover the well known results obtained previously in this 
modeljs^]. It is observed that the components in the 
GS energy are different for the ID Compass and Ising 
models in the transverse magnetic fields for both PBC 
and ABC. It should be pointed out that although the GS 
in PBC and ABC are slightly different in the finite size 
system, they are identical in the thermodynamic limit 
and the essential features in finite-size are also not altered 
qualitatively. Without loss of generality, we will take 
ABC in the following discussion. 



III. FINITE-SIZE SCALING ANALYSIS OF 
FIDELITY AND ENTANGLEMENT 

The GS fidelity and entanglement emerged from quan- 
tum information science have been used in signaling the 
QPTsp, M, E M M M M,M,M- We perform 
finite-size scaling analysis of these two quantities to study 
the criticality of the present model. By using the exact 
GS wave function obtained in Eq.(13), the GS fidelity is 
given by 



Fia,Sa) = |(^o(a)|V'o(a + Sa))\, 



(20) 



where Sa is a small quantity {Sa = 10 ^ is taken in the 
present calculation). Its susceptibility can be written as 



-2lnF 
Xf = hm — - „ . 



(21) 



1,2,. 



N' 



1). Therefore, we 



We calculate the GS fidelity in the (a, /i)-plane, and 
the FS as a function of the transverse field /i for a = 1.0. 
The numerical results are presented in Figs . 3 and 4. 
The absence of the sudden drop to zero of the fidelity 
excludes the level-crossing around the critical point (a = 
±1.0, h=0). In the Kosterlitz-Thouless phase transition, 
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FIG. 4: FS versus ft for J = 1.0 and a = 1.0 with ABC. 



no singularity occurs at the critical point [l8|, so a second- 
order QPT is highly suggested when driving the magnetic 
field, which will be confirmed in the following finite size- 
scaling analysis. 

Next, we illustrate the scaling behavior of average FS 
Xf/N'. The finite-size scaling ansatz for the average FS 
to analyze the second-order QPT takes the form [271 



Xf 



XF 



XF 



^f[N'\h~hr, 



where u is the critical exponent of the correlation length 
and f(x) is the scaling function. This function should 
be universal for large N' in the second-order QPT. As 
exhibited in Fig. 4, the FS reaches a maximum point 
at a certain position hmax- It can be observed in Fig. 
5 that the rescaled FS for larger system sizes tends to 
collapse onto one single curve if adjusting the critical 
exponent v = 1.00 ± 0.02. The scaled average FS at the 
maximum point as a function of N' in log-log scale are 
presented in the inset of Fig. 5. A power law behavior 
-^max ^ j^j-ifj- jg observed in the large N' regime and the 
finite-size exponent extracted from the curve is /i = 2. 
Both values of exponents v and fi in the present model 
are the same as those obtained in the ID transverse-field 
Ising mo del p^. 

We then turn to the quantum entanglement of this 
system. Recently, the concept of concurrence is usu- 
ally adopted as the measure of the local entanglement 
in spin — i systems. The definition of concurrence is 
given by C{i,j) = max[ri{i,j) - r2{i,j) - r3{i,j) - 
r4{i,j), 0], where ra{i,j) are the square roots of the eigen- 
values of the product matrix R = Pijplj in descending 
order (30I [Slj . The spin fiipped matrix pij is defined as 
plj = {cr^ ®(j''^)p*j{cry The is the density matrix 

for a pair of qubits from a multi-qubit state, and has the 



FIG. 5: Finite-size scaling of the average FS according to Eq. 
(22) for hmax — 0, J = 1.0 and a = 1.0 for various system 
sizes N' = 100, 200, 300, 400. The inset exhibits the scaling of 
the maximum of the average FS as the function of log2{N') 
at the critical point (a = 1.0,h=0). 



following form 



3 

a, 0=0 



(22) 

The coefficients are determined by the relations 



Pat3 



tr{(T^<JjPij) 



(23) 



(24) 



According to the refiection symmetry and the global 
phase flip symmetry, considering the Hamiltonian be- 
ing real, the only nonzero coefficients in Eq. (24) are 
Poo,Pii,P22,P33,Po3,P30- Becausc the density matrix 
must have trace unity, so poo = 1- The numerical re- 
sults for the concurrence as a function of h for various 
coupling coefficient a are shown in Fig. 6. It is evident 
that the concurrence gradually increases as enhancing 
a. The minimum of concurrence and a cusp of the first 
derivative of the concurrence occurs right at the critical 
point {a = 1.0, h = 0). 

Furthermore, we calculate the block-block entangle- 
ment both near and at the quantum critical point [2lll38l. 
IsqI . |40| to show the connection with entropy of vacuum in 
the classical conformal field theory in the present model. 
The GS in our model can be completely characterized 
by the expectation values of the two-point correlations 
(aj„a„) — fm,m where m or n is pseudospin site (e.g. 
as,n a2(„_i)+s). Any other expectation value can be 
expressed through Wick's theorem. By eliminating the 
rows and columns in matrix F = fmn, {m,n = 1, 2...N), 
which are corresponding to pseudospins that do not be- 
long to the block, the correlation matrix of the state 
PL is obtained. The corresponding von Neumann entropy 
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FIG. 6: (Color online) The concurrence Ci,2 (left) and the 
derivative daCi,2 (right) versus h for a — 0.6, 1.0, 1.4, J = 1.0 
and N' = 256. 

then takes the form 

L 

Sl = - ^n)log2{l " A„) - Xnlog2Xn], (25) 

ri=l 

where A„ is the nth eigenvalue of the correlation matrix 
Fl- The numerical results for Sl as a function of the 
block size L are presented in Fig. 7. A logarithmic di- 
vergence of 5*^ at the quantum critical point is observed, 
while noncritical entanglement is characterized by a sat- 
uration of Sl for larger L. The coefficient is connected to 
the central charge of the classical conformal field theories, 

limL-^ooSL ~ -^log2L, (26) 

where c = 1 in the compass model. The value of c is dif- 
ferent from that in ID transverse Ising chain (c = 0.5), 
but the same as in ID XX chain without magnetic 
field [111 . It may follow that the block entanglement de- 
pends crucially on the detailed topological structure of a 
system. 

IV. PSEUDO-SPIN CORRELATION 
FUNCTIONS AND MAGNETIZATION 

To explore the essential properties of QPT, we will 
calculate two GS pseudo-spin correlations (erf i erf 2) ^^'^ 
(cTj' j^(t| j^). The numerical results for these two correla- 
tion functions versus a for different magnetic fields are 
presented in Fig. 8. We observe that (cf icrf is a odd 
function of a, while (cf j^crifj) an even one of a. The 
crossing points of (erf 1(7^2) (^1 1^2 1) curves deviate 
a = 1 in the presence of transverse field. The numerical 
results indicate that (trf icrf 2) sensitive to the external 



FIG. 7: The block-block entanglement Sl versus L for J = 
1.0 in the quantum critical point (a = 1.0, h — 0) and the 
noncritical region a = 1.2 with ABC. The inset displays a 
logarithmic divergence for large L at the critical point. 
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FIG. 8: The correlation function with ABC. The parame- 
ters are taken as J = 1.0, A' = 256 and h — 0.1,0.3,0.5, 
respectively. 



magnetic field in the range of a G [— 1 , 1] , but insensitive 
in the other regions. 

As done in Ref. we also calculate the distance de- 
pendence of the pseudo-spin correlator (erf icrf 2+r) with 
ABC for different size in terms of the Hamiltonian (1). 
As shown in Fig. 9 that the correlators at /? ^ 1 de- 
cay in an algebraic way in large r regime, indicating a 
divergent correlation length when approaching the criti- 
cal points. A power law behavior (cf icrf 2+r) °^ 
obtained with 77 = 1.00 ± 0.03, indicating a second-order 
QPT. 

Finally, we calculate the pseudospin magnetization 
(a^) = (erf „) + (cr|„) and the magnetic susceptibility 
X = d{cr^)/dh. The magnetization (cr^) as a function of 
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FIG. 9: (Color online) Distance dependence of {o'2,i<^i,2+r) 
spin correlator at /3 — > 1 for different system size. The pa- 
rameters are J = 1.0, h = and a = 1.0. 
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ID transverse field Ising model, which is given in the in- 
set of Fig. 10 as well. A excellent agreement for the slop 
in the critical regime is clearly shown. 

V. SUMMARY AND DISCUSSION 

By using the method of mapping into a case with 
plural spin sites, we obtain the exact GS energy and 
the GS wave function of ID compass model in a trans- 
verse magnetic field. The pseudo-spin liquid disordered 
ground state is the universal features in the ID compass 
model. Meanwhile, we observe the second-order QPTs 
at (a = ±1, h = 0). The energy gap A will survive even 
in the thermodynamic limit for ft, ^ 0. It is useful for 
supporting naturally robust quantum information. The 
fidelity, the FS, the concurrence, and the block-block en- 
tanglement entropy are also calculated in terms of the 
obtained exact GS wave functions. The finite-size scal- 
ing analysis suggests the second-order QPT occurs by 
driving the transverse field. The pseudo-spin correlation 
functions, the distance dependence of the pseudo-spin 
correlators and the magnetization are also calculated. It 
is observed that the distance dependence of (cf icrf 2+r) 
correlator displays a divergent correlation length when 
approaching the critical points. The obtained scaling ex- 
ponents are nearly the same as those in the ID trans- 
verse field Ising model, suggesting that these two models 
share the same universality class. The scaling exponent 
c = 1.0 of the block-block entanglement entropy is the 
same as the critical XX chain with no magnetic field, 
which shows the different topological structure from the 
quantum Ising model. For the 2D compass model with a 
transverse field, the degeneracy of GS is removed because 
of the destruction of the symmetries. It is expected that 
the QPT becomes more weak and the first-order QPT is 
unlikely. 



FIG. 10: (Color online) The pseudospin magnetization 
((cr^) = (cf^n) + (o"2,n)) versus the transverse field for a — 1.0. 
Inset gives the scaling of the value of IxC^) ~xihc)\ as a func- 
tion of I ft — hc\. The red line denotes that in ID transverse 
field Ising model. 

the transverse field h for a = 1.0 and N' — 256 is exhib- 
ited in Fig. 10. The magnetic susceptibility |x(ft) — x(/ic) I 
versus ftd shows a power law behavior. The exponent 
7 is estimated to be 1.78 ± 0.05 by the slop. It is inter- 
esting that it is very close to the magnetic susceptibility 
exponent 1.75 in 2D classical Ising model. According to 
Eq. (17), we can also plot the similar scaling curve for 
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